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Atmospheric density fluctuations induce phase 
perturbations that degrade the spatial coherence of a laser 
beam propagating through the atmosphere. These 
degradations spread the laser beam and alter the centroid 
and intensity profile stochastically. Turbulent conditions 
are found at virtually all levels of the atmosphere. A 
substantial fraction of the optical turbulence along a 
vertical path arises from the heat flux between the 
atmosphere and the Earth^s surface. This type of 
turbulence is typically within the first 100 to 200 meters 
above the surface. 

During this thesis research, a high frequency acoustic 
sounder was developed to analyze this turbulent layer. The 
primary focus was the development of the command and 
control software required to coordinate the data collection 
and reduction. The system was used at two sites and should 
prove useful in quantifying the effects of optical 
turbulence within the surface boundary layer on laser and 
optical system performance. 
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I. INTRODUCTION 



A coherent laser beam propagating through the 
atmosphere is very susceptible to numerous turbulence 
inflicted degradations. As electromagnetic waves transit 
the turbulent regions, atmospheric irregularities randomize 
the amplitude and phase of the wave. In order to quantify 
the altitude dependence of the atmospheric turbulence, a 
high resolution turbulence profiler is needed. Acoustic 
echosounders are frequently used to detect and measure 
atmospheric density and velocity irregularities resulting 
from air currents, temperature inversions, humidity 
variations, mechanical turbulence and other density 
fluctuations . 

Presently two atmospheric optical parameters, the 
spatial coherence length (ro) and the isoplanatic angle 
(00) are measures of the perturbation of an electromagnetic 
wave propagating through the atmosphere, and are accurately 
measured by optical systems developed by Walters [Refs. 1 
and 2]. Although these systems measure a path integral of 
atmospheric turbulence with high accuracy, a major drawback 
of each system is that no provision is available for 
determining the height of the atmospheric disruptions. If 
these disturbances are found to exist very near the 
surface, it may be possible to negate their effect by 
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elevating the optical systems or controlling the generation 
of turbulence. 

This thesis deals with the design, construction and 
implementation of a high frequency, acoustic echosounder 
which will accurately analyze the atmospheric density 
fluctuations within approximately 200 meters of the 
surface. As this project is a product of the research and 
efforts of two students, the work was appropriately 
divided. My particular task was to devise a computer 
program which would control all aspects of the echosounder 
operation. Areas of particular interest were the 
integration of hardware and software within the data 
acquistion system, developing plotting algorithms, 
controlling the system timing, and setting the input 
parameters which determine the range and sensitivity of the 
device. My colleage, LT . Weingartner [Ref. 3], dedicated 
his efforts toward the actual design and hardware 
development of the echosounder. 

Acoustic echosounders have been developed and in use 
for many years and have proven to be valuable probes for 
analyzing the structure and dynamics of the lower 
atmosphere [Refs. 4 through 7]. Devices similar to ours 
have been used to obtain profiles of the atmospheric 
density and temperature fluctuations [Refs. 8 and 9]. Our 
device used a high speed HP217 computer to control and 
monitor the echosounder which provided real time 
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information. Additionally, we had the ability to store and 
reproduce the atmospheric profile plots at will. This 
information was crucial in analyzing the atmospheric 
turbulence degradations on laser and electro-optical 
system’s performance. 

This document shall address the theoretical background 
of acoustic sounder operation in Chapter II. Chapter III 
shall present a synopsis of the system design and 
associated software. A summary of actual data collected at 
two sites is analyzed in Chapter IV and Chapter V discusses 
the conclusions and recommendations found as a result of 
our research. 
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II. BACKGROUND 



Acoustic echosounders probe the atmosphere by 
transmitting a pulse of acoustic power which is 

subsequently scattered back from the atmosphere by 

temperature and velocity inhomogeneities. The echosounder 
(echosonde) equation, often referred to as the radar 
equation in meteorology, is used to determine the 

backscattered acoustic power. This equation is summarized 
by Neff in Reference 10 and is based upon the work of 
Tatarski [Ref. 11] and Little [Ref. 5], 

Pr = Er[ PtEt ][ exp(-2aR) ][ ao(R.f) ][ J^cr ][ AGR- 2 ] 

where 

Pr is the electrical power returned from a range R, 

Pt is the electrical power transmitted at frequency f, 
Er is the efficiency of conversion from acoustic power 
to electrical power by the transducer, ' 

Et is the efficiency of conversion from electrical 
power to acoustic power by the transducer, 

exp(-2aR) is the round trip power loss due to 
attenuation where a is the average attenuation (meters' i ) 
to the scattering volume at the range R (meters), 

ao(R,f) is the scattering cross section per unit volume 
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at a distance R and frequency f, 

c is the local speed of sound (meters/second), 

X is the pulse length (seconds), 

A is the aperture area of the antenna (meters^ ) , 

R is the range (meters), and 

G is the effective-aperture factor of the antenna. 
Empirically measuring or calculating the values for all 
other terms, one can use this equation to determine 
0 O(R,f), the scattering cross section per unit volume; that 
is, the fraction of incident power backscattered per unit 
distance into a unit solid angle at a frequency f. Based 
upon experimental results, Tatarski [Ref. 11] expresses the 
acoustic backscatter cross section per unit volume, 
ao(R,f), in the equation. 



Ct ^ 

cro (R,f ) = 0.0039 kW3 

To = 



where 

k = 2x/ A is the incident acoustic wavenumber at 

wavelength A, 

To is the local mean temperature in degrees Kelvin, and 
Ct^ is the temperature structure parameter. 

Combining this equation with the echosonde equation, one 
obtains a volume-averaged measure of Ct®. 
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11 2 1 Pr 

Ct^ = To^ k-i/3 R^exp(2aR) 

0.0039 ErEt ct AG Pt 

Hall and Wescott [Ref. 12] calculated a beam-shape 
compensation factor of 0.40 for a piston source antenna 
with a uniformly illuminated square aperture. This value 
is the same as the effective-aperture factor, G, and can be 
substituted into the above equation. Approximating the 
aperture area of the antenna to be equivalent to 25 times 
the aperture area of a single speaker having a diameter of 
7.620 centimeters, we get a value of 0.1140 square meters. 
Combining these values with the numerical constants in the 
above equation, we can simplify the equation for Ct^. 



1 1 Pr 

= 11245 To= k-W3 R^exp(2aR) 

Er Et ct Pt 



The acousto-electrical efficiency factors for our 
echosounder were measured with a calibrated microphone and 
an anechoic chamber. The calculated values for Er and Et 
are each 0.5. Substituting these values into the above 
equation, we can further simplify the equation for Ct^ . 



To^ Pr 

Ct= - 44980 k-i/3 R2 exp(2«R) 

CT Pt 
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The above equation was incorporated into the computer 
program that provides the system control, data acquistion 
and data processing techniques [Appendix A] . The- reduced 
data was then used to measure the temperature structure 
parameter as a function of time of day and altitude for 
various sites. 
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III. SYSTEM DESIGN AND EQUIPMENT DEVELOPMENT 



A. HARDWARE 

Most of the hardware used in this system was standard 
scientific equipment and is illustrated in Figure 1. The 




Fig. 1. Echosounder Equipment Arrangement 
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pieces which were specially designed by LT Weingartner 
|Ref. 3] are the acoustic array and the enclosure for this 
array. The design considerations for the speakers' and the 
array format are given in Appendix B and the design -'»f Lhe 
enclosi.ire to house this array is outlined in Appendix (h 
Aside from these two pieces of equipment, there are 
primarily six other components which comprise this 
echosounder. Below is a brief description of each of these 
additional components . 

^ • HP-2 0 0_^ S e ries Co m p uter 

The Hewlett-Packard (HP) 20U Series Computer 
includes a 20 megabyte hard disk drive, a floppy disk drive 
and an associated printer and monitor. The HP200 Series 
Computer was the central control component for the entlr-^ 
echosounder arrangement. The computer used was an HP21/ 
pr‘'^grammed in Basic 3.0 and equipped with an Infotek 
Basic Compiler and an Infotek FP210 Floating P.>int 
A'-^ce I era tor to enhance the speed of pi-ogram execution. The 
prr^gram defined all parameters for the HP3314A Function 
Generator as well as execute the trigger command which 
produced the echosounder transmitted acoustic pulse. The 
computer also received the data from the acoustic array via 
the pre-amplifier, bandpass filter and analog to digital 
converter. The computer then conducted the data redu-^tion 
routines to produce and display a high resoiuti''n 
atmospheric profile . 
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2 . HP3314A Function Generator 



The Hewlett-Packard (HP) 3314A Function Generator 
was a multimode function generator capable of providing 
sine, square and triangular wave functions as well as any 
desired waveform ranging in frequency from 0.001 Hertz to 
19.999 Megahertz. The HP3314A Function Generator was used 
to supply a pulse of an integer number (usually 100) of 
sinusoidal cycles of constant amplitude to the QSC Model 
1700 Audio Amplifier. A constant frequency setting of 5000 
Hertz was used for all data runs. 

3 . QSC Model 1700 Audio Amplifier 

The QSC Model 1700 Audio Amplifier is a high power 
amplifier which can supply 350 watts over an 8 ohm load. 
This amplifier was used to boost the output voltage of the 
function generator by a factor of 20 from 1 . 5 volts to 30 
volts. The power supplied to each speaker in the array 
during operation was then 37.5 watts. 

4 . Pre-amplif ier 

The pre-amplifier was designed and constructed by 
Walters [Ref. 13] to supply a gain of 1000 to the returned 
signal and isolate the data acquistion components from the 
transmitted pulse. The pre-amplifier may be thought of as 
a safety and switching mechanism for the system. An LT 
1037 and an LT 1007 Operational Amplifier were selected for 
use in the pre-amplifier based upon their low noise 
properties which were evaluated by physically incorporating 
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them into the device and measuring their noise and gain 
characteristics with an HP3561 Spectrum Analyzer. The 
diodes used are IN 4000 Series Rectifiers capable of 
carrying one ampere. These rectifiers serve to isolate the 
power amplifier from the amplifier input stage and to limit 
the voltage applied to the operational amplifiers. The 
electrical diagram of this device is outlined Figure 2. 




Fig. 2. Schematic of Pre-amplifier 
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5 . Rockland Wavetek Model 852 Filter 



The Rockland Wavetek Model 852 Filter operated as a 
48 db per octave bandpass filter to suppress the broadband 
noise of the system. High and low bandpass settings of 
5500 Hertz and 4500 Hertz respectively were used. The 
filter response at these settings is illustrated in Figure 
3. 



START I 000 Hz range: -5 dBV STATUS: PAUSED 

A:NAG RMS:50 




START: i 000 Hz BU: 477.43 Hz STOP: 51 000 Hz 

X: 3000 Hz Y:-20.24 dBV 



Fig. 3. Bandpass Filter Response 

The three db (half power) bandwidth is approximately 1 KHz 
wide. Although it would be desirable to reduce the 
bandpass to between 50 and 100 Hertz without impairing 
signal response, the Rockland Wavetek Model 852 Filter was 
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"the best filter available. It performed adequately in 
spite of the large IKHz bandwidth. 

6 . Infotek AD2QQ Analog to Digital Converter 

An Infotek AD200 12 bit Analog to Digital Converter 
was used to digitize the signal voltage for the computer 
system at a 12.5 to 20 KHz sample rate. 

B. SOFTWARE 

The software and HP200 Series Computer system were 
responsible for controlling and monitoring every 
operational phase of the hardware components. The program 
which accomplished this task was entitled "ACRDR" and is 
listed in Appendix A. This program was modeled after a 
program written by Walters [Ref. 13] but each program 
performs a distinctly different computational task. This 
program was written in HP Basic 3.0 and was compiled by an 
Infotek BC203 Basic Compiler to enhance the speed of 
execution . 

The program "ACRDR" is easily broken down into a number 
of blocks and subroutines which performed specific 
operations. These sections are outlined in a flowchart 
(Fig. 4). Each block is straightforward in its purpose, 
and the program is designed to be as helpful to the user as 
possible. As a prologue to the actual program code, there 
is a listing of all the program variables with a short 
description of their use. Such a listing familiarizes a 
user with the computations to be made and also provides 
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quick references 



for any 



future modifications 



to be 



implemented. 




Fig. 4. Flowchart of Fundamental Program Operations 

The initialization procedures which set all the 
parameters necessary for data collection follow the 
variable definition section. These parameters control such 
things as the contrast between background and return signal 
and the setting of the computer’s internal clock. All 
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options of the HP3314A Function Generator are decided by 
the user at this phase. Another initialization involves 
preparing the Infotek AD200 Analog to Digital Converter for 
operation. Finally, internal arrays are dimensioned, 
integer variables are defined and the computer function 
keys are redefined to suit the "ACRDR" program. 

Now that the system is ready to operate, certain 
quantities must be determined based upon atmospheric 
conditions at the time of data collection. Fi^om the 
ambient temperature, a speed of sound in air is determined 
by the equation in Kinsler, et al , [Ref. 14: p. 106] 

C ^ Co [ 1 + (T/273) ]i/2 , 



where 

C is the local speed of sound, 

Co is the speed of sound in air at O'" Celsius (331.6 
m/s ) , and 

T is the local temperature in degrees Celsius. 

Using this speed with a user decided sampling frequency and 
number of data points, the maximum range of the traces is 
determined by the relation, 

RANGE = C/2 • l/( SAMPLE FREQUENCY) • N , 
where 

RANGE is the range of the echosounder in meters, 

SAMPLE FREQUENCY is the user desired rate of data 
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acquisition, and 

N is the number of data points per trace. 

Typically, the vailues of many of the parameters are not 
varied. The value of N above is almost always taken to be 
the allowed maximum of 16301 data points and the speed of 
sound can be safely estimated to be 340 m/s at room 
temperature. By using a typical sample frequency of 20000 
Hz, a range of approximately 135 meters is obtained. 

Another parameter determined is the atmospheric 
attenuation coefficient. This calculation is done in a 
subroutine obtained from Reference 15 and is necessary in 
this program for the calculation of Ct 2 . Other 
computations made prior to program execution include the 
pulse length, wavenumber, and the D-C offset of the 
equipment . 

At this point, the computer is finally ready for data 
collection. The screen setup displays a distance versus 
time plot and the internal clock of the computer is 
synchronized with the pulse from the HP3314A Function 
Generator. Following transmission of the pulse, the return 
signal is received, digitized, stored in an array and data 
reduction commences . A block averaging technique is used 
in which a block of data points the size of the number of 
emitted cycles is summed and averaged. Returns are plotted 
as darkened areas on the echosounder traces with the 
intensity of the darkened area linearly related to the 
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magnitude of the return signal. This same procedure is 
used in the computation of Ct 2 with some very important 
differences. Aside from simply correcting for Tir^ 
spherical divergence, Ct 2 is corrected for electronic gain, 
the ratio of power returned to power transmitted, the 
efficiency of the speakers as transmitters and receivers, 
the area of the speaker array, the atmospheric attenuation, 
the pulse length, the temperature, the scattering cross 
section per unit volume at a specific range and frequency 
and finally the effective aperture factor of the antenna. 
Addit ionally , Ct 2 (Range) is averaged for a particular 
altitude over 15 minute intervals. Finally, after each 
pulse is reduced, a corresponding mean square noise level 
is determined. This is done by averaging each block 
average at maximum range until at least ten values have 
been used in the average. An upper limit of five over the 
average noise figure (this corresponds to voltage 
fluctuations on the order of 10"'^ volts) is set on the 
routine to avoid averaging any strong return signals or 
anomalies such as passing aircraft. After the ten values 
are averaged, every subsequent pulse is averaged into all 
the preceding noise levels and removed from each subsequent 
return signal. 

Ultimately, after each 15 minute interval, or at the 
users request, the display terminal image is printed. At 
this time the Ct 2 computations are conducted and plotted on 
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the screen and printer as range versus logi o of Ct 2. Upon 
completion of the printing of these plots, the system 
begins the data processing for the next 15 minute interval. 

There are certain options built into the program to 
allow a user to change various aspects of operation. The 
function keys allow the user to change the local 
temperature, sample frequency or intensity factor during 
program execution. Additionally, the user can quit or 
restart the program, print the partial trace on the screen 
or elect to save a trace on a floppy disc. The save 
routine is invoked for the subsequent 15 minute interval 
after the appropriate function key is depressed. Saving a 
future trace may seem awkward, especially if a user would 
like to keep an interesting trace which is presently on the 
terminal. This problem cannot be readily solved unless 
each trace is recorded to disc without user intervention. 
At present, this is not done because only 8 traces (2 hours 
of data) can be written to a floppy disc before it is full. 
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IV. 



A. ECHOSOUNDER PERFORMANCE 

An analysis of the acoustic echosounder output was 
conducted to determine the validity of previously 
determined echosounder parameters. The e~ i decay time 
constant for applied voltage was calculated in Appendix B 
to be approximately 900 ps or for convenience, 1 ms. 
Reviewing typical echosounder traces and Ct 2 plots 
[Appendix D] , it was evident that the recovery time was 
typically found to be on the order of 33 ms or roughly six 
meters past the end of the transmitted pulse length of 3.4 
meters. The recovery time is consistent with the time 
required for the 30 volts on the drivers to decay to the 
microvolt level. A more detailed analysis of the hardware 
is found in Reference 3. 

B. SITE EVALUATION 

Echosounder data was collected at two different 
locations. The primary data collection site was the upper 
roof of Spanagel Hall at the Naval Postgraduate School, 
Monterey, California. This site was chosen simply for 
convenience. Data gathered at this location is believed to 
represent the California coast during the spring near sea 
level. The second site chosen was in the vicinity of the 
24 inch telescope at Lick Observatory, San Jose, 
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California. This site is located atop Mt. Hamilton at an 
altitude of approximately 5700 feet and nearly 20 miles 
inland from the coast. 

These two data collection sites represent areas of 
differing atmospheric air pressures, water vapor pressures, 
local temperature ranges, and local wind velocity ranges. 
These characteristics all play important roles in effecting 
the local atmospheric turbulent conditions and thereby the 
atmospheric structure parameter, Ct2. 

In addition to collecting echosounder data at Lick 
Observatory, simultaneous measurements of the isoplanatic 
angle (00 ) and spatial coherence length (ro) were made with 
systems developed by Walters [Refs. 1 and 2]. A basic 
knowledge of these two systems is necessary to understand 
the correlation procedures made. The isoplanatic angle 
(00) is primarily an upper atmospheric measurement which 
indicates atmospheric disruptions at a range of 2 to 15 
kilometers. The spatial coherence length ( ro ) is a measure 
of the effects of the entire atmospheric blanket on 
coherent light transmission. A close comparison of all 
three data sets should give us an accurate description of 
both the lower and upper troposphere as well as the 
stratosphere above. 

Based upon isoplanatic angle (0o) and spatial coherence 
length ( ro ) measurements at Mt . Wilson in California, a 
strong correlation between the two measurements occurs if 
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the low altitude boundary layer contribution is 
sufficiently small. This strong correlation helps to 
reinforce the overall description of the atmosphere at the 
time of data collection. Figures 5 and 6 graphically 
illustrate the atmospheric measurements made at Mt. Wilson 
on 2 April 1987 by Walters. 




Fig. 5. Isoplanatic Angle Measurements, Mt. Wilson 
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500 



2 Rpril 1387 Mt Wilson 




Fig. 6, Spatial Coherence Length Measurements, Mt. Wilson 

The strong correlation between the isoplanatic angle 
( 00 ) and the spatial coherence length (ro) is especially 
evident between the hours of 0700 and 1300 universal time. 
The close tracking of these two measurements during this 
time interval indicate that the upper atmospheric 
conditions, as measured by 9o , are dominating the entire 
atmospheric profile as measured by ro . Unfortunately, the 
existence or non-existence of any turbulent surface effects 
cannot be ascertained by the employment of the above two 
systems alone. However, employment of these two systems 
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together with the acoustic echosounder should enable us to 
produce a complete atmospheric profile with strong 
correlation between all three atmospheric measurements. 

On 9 and 10 April 1987, all three systems were operated 
at Lick Observatory. Again, a good correlation between the 
isoplanatic angle (0o) and the spatial coherence length 
( ro ) measurements was noted. In addition, a strong 
correlation between the spatial coherence length ( ro ) and 
echosounder measurements was present. A comparison of the 
atmospheric data in Figures 7 and 8 shows a good 
correlation between the two parameters especially during 
the 09.30 to 1130 time interval on 9 April 1987. However, 
during subsequent hours the isoplanatic angle (9o) 
measurements remain high (=12 prad) indicating relatively 
calm turbulent conditions in the upper atmosphere while the 
spatial coherence length (ro) values drop sharply aftei' 
12:00 Universal Time indicating dominant and increasing 
lower atmospheric turbulence. This increase in the low 
level turbulence should be evident in the echosounder data 
commencing around 1200 universal time (0400 local standard 
time) on 9 April 1987. A comparison of the echosounder 
data in Figures 9 through 11 illustrates this increase in 
the local surface turbulence. 
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Isoplanatic angle (urad) 



20 



3 Rpr i 1 1937 Lick 




Fig. 7. Isoplanatic Angle Measurements, Lick Observatory 




Fig. 8. Spatial Coherence Length Measurements 
Lick Observatory 
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9 Bpr 193 7 LICK OBSERVFITORY 




TIME LOCRL (STD) 

Fig. 9. Echosounder Trace, Lick Observatory 



9 Rpr 1987 LICK OBSERVATORY 
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4j00 4:05 4:10 4:15 



TIME LOCAL (STD) 

Fig. 10. Echosounder Trace, Lick Observatory 



31 



9 Apr 1987 LICK OBSERVATORY 

135 r .... ..... 




5:00 5:05 5:10 5:15 

TIME LOCAL (STD) 

Fig. 11. Echosounder Trace, Lick Observatory 

The strong correlation between the echosounder data around 
4:00 to 5:00 Standard Time and the ro measurements around 
12:00 Universal Time combined with the lack of correlation 
between the 9o and ro measurements indicate that the lower 
atmospheric and surface turbulence are dominating the 
atmospheric profile during this time period. 

Data collected on 10 April 1987, again illustrate the 
strong correlation between the three atmospheric 
measurements made. A comparison of the 0o and ro 
measurements during the time interval of 0700 and 1300 
Universal Time indicate steady turbulent conditions in the 
upper atmosphere and greatly varying turbulent conditions 
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at lower atmospheric levels. This is evident in Figures 12 
and 13 by the consistent values of 00 during the time 
period compared with the steady increase and eventual 
decline of the ro values during the same time interval. 
The trace variations in Figure 13 during the hours of 0800 
and 1200 are indicative of a period of decreasing lower 
atmospheric or surface turbulence followed by the onset of 
an increasingly turbulent period around 1200. This 
turbulent trend is strongly supported by the echosounder 
data in Figures 14 through 21. 




Fig. 12. Isoplanatic Angle Measurements, Lick Observatory 
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RRNGE IN METERS 




Fig. 13. Spatial Coherence Length Measurements, 

Lick Observatory 



10 Apr 1387 LICK OBSERVATORY 




TIME LOCAL (STD) 



Fig. 14. Echosounder Trace, Lick Observatory 
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RHNGE IN METERS RRNGE IN METERS 
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Fig. 15. Echosounder Trace, Lick Observatory 




TIME LOCAL (STD) 



Fig. 16. Echosounder Trace, Lick Observatory 
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RANGE IN METERS RANGE IN METERS 



10 Apr 1987 LICK OBSERVATORY 




TIME LOCAL (STD) 



Fig. 17. Echosounder Trace, Lick Observatory- 



10 Apr 1987 LICK OBSERVATORY 




3:00 3:05 3:10 3:15 

TIME LOCAL (STD) 

Fig. 18. Echc ounder Trace, Lick Observatory 
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RRNGE IN METERS RRtvJGE IN METERS 




TIME LOCAL (STD) 

Fig. 19. Echosounder Trace, Lick Observatory 
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Fig. 20. Echosounder Trace, Lick Observatory 
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10 Apr 1987 LICK OBSERVATORY 




TIME LOCAL (STD) 

Fig. 21. Echosounder Trace, Lick Observatory 

The large 200mm coherence lengths around 11:00 Universal 
Time are consistent with the low turbulence evident in the 
echosounder profiles around 3:00 Standard Time followed by 
a pre-dawn increase in the surface turbulence. 

Only acoustic sounder atmospheric measurements were 
made on the upper roof of Spanagel Hall at the Naval 
Postgraduate School. Data runs on 26 and 27 April 1987, 
included both echosounder data and the associated Ct 2 
plots. A representative sample of the data collected 
during this period is included in Appendix D. Of 
particular interest is the data in Figure 22, which shows 
the maritime boundary inversion layer at about 100m and 
being perturbed by convective plumes at lower altitudes. 
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Fig. 22. Echosounder Trace and Ct 2 Plot, NPS Roof 
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V. CQNCL.USJQNS. ■ AND-RE.C£>MMENDAT1QNS 



A high frequency ( 5KHz ) phased array echosounder 
constructed to measure low level turbulence appears to work 
within the 10 to 135 meter altitude range. When using the 
device, detailed profiles of the short range atmospheric 
density fluctuations were obtained. The profiles were 
found to correlate very well with the measurements of the 
isoplanatic angle (0o) and the spatial coherence length 
( ro ) during periods of simultaneous operation. This short 
range echosounder, when used in conjunction with the other 
atmospheric measuring devices, is an invaluable tool. The 
graphical output provides us with a more complete 
description of the atmosphere that can be used to calculate 
Cn® , the atmospheric index of refraction structure 
parameter . 

Further software development such as the incorporation 
of Fast Fourier Transform routines into the echosounder 
program will determine the radial velocity profile of the 
return signal. Such an improvement will determine the 
velocity of the probed air masses and plot these 
echosounder traces as a function of color intensity. 

Other areas of further reasearch include the design and 
testing of different array patterns. One such pattern, the 
hexagonal array, is already undergoing tests. As the 
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speaker array evolves, it is also evident that the 
enclosure must follow suit to accomodate the new array 
design and associated beam patterns. Further improvements 
in software routines are also inevitable. It is always 
desirable to store on a disc all data collected at a site. 
Presently, the floppy disc capacity is inadequate for the 
storage of more than two hours of data. Furthermore, if 
FFT and Doppler routines are to be added, the computational 
speed of the computer will inhibit the pulse repetition 
rate. A special data acquisition technique called "Direct 
Memory Access" may then have to be introduced into the HP 
computer to allow simultaneous data collection and 
processing . 

Finally, this product may be used in various 
applications not already addressed. One such use may be to 
measure windshear at airports . Other researchers have 
expressed interest in using the echosounder to measure 
arctic atmospheric conditions during meteorological 
surveys. Assuredly, as the device evolves into a more 
compact and highly efficient instrument, its range of 
application will continue to expand. 
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APPENDIX A 



ACOUSTIC ECHOSOUNDER PROGRAM 



10 
20 
30 
40 
S0 
60 
70 
80 
90 
100 
1 10 
120 
130 
140 
1S0 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
L. i XJ 
288 

280 

Tn>n> 

U'U'U 

318 

•7''^ n, 
U tJ 



350 



410 






» ACRDR ♦ 

* COMPUTER SOFTWARE WRITTEN BT LT. M. WROBLEWSKI FOR AN » 

» ECHOSOUNDER BUILT FOR A MASTERS OF SCIENCE THESIS • 

» BY LT. M. WROBLEWSKI AND LT. F. WEINGARTNER » 

♦ ADUISDR: PROF. D. WALTERS * 

» JUN. 1987 » 






! The conputer prograi^ receives infornation fron an acoustic 
! array through an A“D converter. This information is the 
! returned signal of an acoustic pulse as it passes through 
! the atmosphere. The data is then used to display the return 
! intensity with distance as a function of time. This is a 
j short range device (from 150m to 200m). 

! 

j 

I LIST OF VARIABLES 
! 
j 

I I S K u J counters for loops 

j 

! Rec_num & Nrec - real and integer represent at ion of the 
! record number for storage 

j Plotnum - counter used to insert form feeds between plots 

! 0 i 5 C__addres s3 — storage location of data file 

} 

1 r i 1 e 3 
! m 3 u Num 



Point ( * ) 



I Da t ( " ) 
! 



...... I 1. .1... . .l.X. 

t ixc Tioi'ic uacu uu slOI e Ltaua 

string and real representat ion of the riumber of 
cycles and the number of points used in the 

computation of the block average 

- real number representing the noise and range 
corrected average over given number of points 

array of a-d converter output after sampling 
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A nm 
^C<J 


! D2( » ) 




array of reduced and averaged data including 


430 


! 




offset, correction for noise, and range 


/< < rn 


! 




correct ions 


450 


! 






4G0 


! Hrs - 




the integer hours 


4 nr\ 

^ { <j 


! 






480 


! Min 




the nunber of ninutes 


490 


! 






500 


! Qtrhr 


2, 


Qtrnin - keeps track of the passage of each 


510 


! 




15 minute interval 


520 


1 






530 


! NumI 


- 


counter for copiputation of block averages 


540 


! 






550 


! Num2 


- 


counter for the number of block averages made 


5G0 


! 






570 


! M - 


the counter used in the plot label routines 


500 


! 






530 


! TiMcnt 




time interval between data samples in 


500 


j 




nanoseconds 


610 


! 






520 


! 






630 


! Knax 


- 


total number of block averages computed 


640 


! 






650 


1 FreqS 


& 


Freq - the frequency input of the HP3314A in kilohert 


660 


j 






670 


! Zones 


- 


the appropriate time zone the operator desires 


600 


! 






690 


! Sites 


- 


the name of the appropriate site of data collection 


700 


j 






710 


! Save_plt 


& Cntrl - on/off toggles used to determine whether 


{ iUXJ 


1 




a particular run will be saved to a disc 


n"7f> 
i yJMJ 


! 






740 


! Npiot 


- 


the nupiber of increment along the entire 


750 


1 




horizontal axis 


7G0 


! 






770 


! Offset 




the computed D -C offset for the system prior 


nnr% 
{ O'a 


! 




to data acquisition 


790 


! 






000 


! Noise 


- 


a running total of the noise accumulated at the 


n 1 r\ 
U 1 KJ 


} 




maAir'iur'i » Siiye? useu ^.o i Imu aiiu rei'’iOve the 


020 


1 




a V e r a g e noise 


030 


1 






840 


! Liput 


- 


the upper bound on the noise figure used to insure 


850 


! 




that large returns are not included in the 


860 


1 




noise computations 


370 


1 






880 


! Ki - 


a 


counter of the number of traces used to determine 


890 


1 


the noise average 


900 


j 
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910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1000 

1090 

1 100 

1110 

1120 

1130 

1 140 

1 150 

1 160 

1 170 

1 180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1200 

1290 

1 300 

1 -7 ^ r\ 

I U i 'U 

A ~rn n 

I o/ 

1330 

] ^4*0 

I3S0 

1360 

1370 

1330 

1390 

1400 



I Sanave 

j 

I 



a sample average used in the computation of the 
D-C offset 



! Off^ave 
j 
! 

! Counts Sf 
! 

! 

! Samfreq 



a toggle used to inhibit the firing of the 
HP3314A during D-C offset computations 

Cnt - the string and real represent at ion of 
the number of data points to be taken 

the user input sampling frequency desired 



Tminc - 1 over the sampling frequency; the time between 

samples 

Times - the string required by the A-D converter to 
sample at the desired rate 

COkelvin - the speed of sound at 0 degrees celcius in m/s 

Temp - temperature of the surroundings in degrees celcius 

Spd - the speed of sound in air computed for the input 
temperature 

Rspd - the relative speed of sound that the echosounder 
sees which is half the computed speed 

Rdist - the distance traveled in one time increment 

Maxrng - the maximum range of the echosounder found 

by multiplying the distance per time increment 
by the number of time increments 

Exrng - the range rounded to the nearest value evenly 
divided by 15; used solely for plotting 
purposes 

Far ' used to estimate the far field; this value keeps the 
range correction from being applied to data very near 
to the source 

1 1 V 1 ^ the user input intensity level divisor; this value 

sets the screen contrast in data return 



Seal e y 



the ratio of exrng to maxrng; this value is used 
to keep the plotted range of data accurate 



Npoint u Npointl & Npoint2 - variables used to keep track of 

time passage between pulses 
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1 410 

1420 

1430 

1440 

1450 

14G0 

1470 

1430 

1 490 

1500 

15T0 

1520 

1530 

1540 

1550 

I5G0 

1570 

1580 

1530 

1G00 

1G10 

1620 

1G30 

1G40 

^G50 

1GG0 

1G70 

1G80 

1G30 

1700 

1710 

i nn r% 

\ ( tUXJ 

i r 

i n A rk 
I ( 

1750 

I7G0 

1 nnm 
] t i ^ 

1 ^nn, 

I ( 0<J 

1790 



< rv 4 
1 C 1 o 

1 on rJk 
\ '-J 

1 n*r/Tk 
i 

i n * n> 
I 

1850 

18G0 

1 nnr\ 
\ ^ t XJ 

1 

I oo-u 

1830 

1300 



! TZ Sr T1 Sf T0 - used to synchronize data collection yith the 
I clock 

! 

I TIMEDATE - the internal clock of the computer 
{ 

! Dvd - divisor of the data average; used because computer 
! multiplication is faster than division 

! 

! R - the range of a block of data samples 

I 

! Run_ave - the running sum of the block samples 
! 

! X *- the value of the data point less the D~C offset 
! 

! Ns - the final average of the block of data points; also 
! used in the computation of the noise figure 

! 

! Corr - the noise correction applied to the data samples 
! 

I Time - the horizontal position of the trace on the plot 
I 

I Timedist - the horizontal width of the trace on the plot 
I 

! Endjt - the last point of the plot vertically taking 
I into account the scaling factor sealery 

I 

! Inc - The vertical increment along the plot 
I 

I Z - The final reduced data points which are output on the plot 
I 

! Cis - The vertical height of the trace on the plot 
! Ap$ ,Frq$ ,Nm$ ,En^ ,VoS ,Hz$ - strings needed to set the HP3314A 
! Amp$ - userv input amplitude for the HP3314A 
1 Pl5_lng - the pulse length of the burst 

! Ct(») - the atmospheric temperature structure parameter 

! Y the number or traces used in computing Ct( ' ) 

! Ypl “ the vertical position on the Ct plot 

! Pnt “ the horizontal position on the Ct plot 
I 

! K3 - the wave number to the 1/3rd power 
I 

! A - The area of the receiver (array) 

! 
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1310 


[ 


G - 


The effective aperature factor 


1920 


1 




' 


1 930 


1 


Er - 


efficiency of coversion of acoustical power 


1940 


! 




to electrical power on the recieve side 


1350 

13B0 


! 

! 


UJ 


the efficiency of conversion of electrical 


1370 


! 




power to acoustical power on the transr-Ut side 


1380 
1 330 


1 

! 


Pt - 


The computed transmitted power to the 


nrsrso, 

L.XJ'O'U 


i 




acoustic array 


ntyin, 
i.'O 1 XJ 

2020 


j 

! 


Gain 


the electronic gain of the equipment 


2030 

2040 


1 

1 


Z imp 


the speaker array impedance 



20S0 

n nn r% 
L. ub xJ 

2070 
2080 
2030 
2! 00 
2110 
2120 
2130 
2U0 
2150 
21B0 
2170 
2180 
2130 
2200 



OPTION BASE 1 
! 

1 

initialize the arrays Sr set dimensions 
declare all integer variables 



DIM Disc_address$C20] ,File$C30] , Point ( 300 ) ,Ct( 200 ) 
INTEGER I ,Hr .Plotnum ,Pr int_key ,Num ,Numl ,Num2 ,M ,K 
INTEGER Rec_num ,Kmax 
INTEGER 02(300), Dat(lB301 ) BUFFER 



20 I 



initialization rout ines .... set time, set HP3314A Function 
Generator 



Rstrt:CALL Freq_init ( FreqS ,N$ ) 



22B0 


CALL Init_ad200 




nnn n, 
C. i. l HJ 


CALL Set_t ime( Zone$ ) 




2230 


INPUT "SITE NAME ",Site$ 




2230 


Gave p 1 1 =0 




^ U'U^ 


0 li k. ( X - 'o 




^ -7 1 rx 

i_ ^ 1 ^ 


1 




2320 


! 




n~7-7rs 

CUU'O 


! keyboard set_up ~ sets labels on the 


computer function 


2340 


! keys 




2350 


1 




23B0 


! 




2370 


OUTPUT KBD; “SCRATCH KEYE” ; ! CLEAR THE KEYS 


2380 


CONTROL 2 ,2; 1 




2330 


ON KEY 1 LABEL "NEW TEMP" GOTO 


Speed 
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hj M rj 



2400 

2410 



L. H-OU 
n * A rs 

2450 



n "»ntr 
i, I U< 



ON KEY 2 LABEL “INTENS. FACTOR" GOTO Lvl 
ON KEY 4 LABEL "PRINT TRACE" GOTO Prt_dt»ip 

ON KEY 5 LABEL "SAOE NXT PLOT" GOTO 5vpt 
ON KEY B LABEL "CHANGE SMPL FRO" GOTO Sfrq 
ON KEY 7 LABEL "RESTART" GOTO Ratrt 
ON KEY 3 LABEL "QUIT" GOTO Quit 



n 4 n/n 


1 


n 4 


! 


n 4 nrk 


1 Set CODatanto 


2430 


j 


2500 


! 


2510 


Di5C_addre55$ 


2520 


Ma>;rec=720 


2530 


Nplot-300 


2540 


Plotnun=l 


2550 


Of f 3et=0 


2560 


GINIT 


2570 


No i 5a=0 


2580 


LiMi t= » 0 


2590 


Y=0 


2600 


Kl-0 


2610 


Sanave=0 


2620 


« 1 1 40 


2630 


G=.4 


2640 


Er=.49S 


2G50 


Et=.49G 


2660 


Gain=73 1 2B *'2 


2670 


Zirop=»i 2 . 1 2 


2680 


1 


2630 


I 



: ,700,1 ,0" 



Input and calculation of terns used by the A-D converter and 
the conputation of range 



2740 

2750 

27GO 

n n^r% 
^ t t tj 



2730 

n nn>rx 



n nn rk 
^ U 



PRINT "Enter the number of data points desired (max 1G301)" 
PRINT "Use increments of the number of cycles then add I" 
PRINT “ex: <150 cycles X 100 data point s/cyc le )i- 1 = I 500 1" 
PRINT "RECOMMEND USING MAX UALUE OF 1B301" 

INPUT Counts 
! 

! uuU|jUl [juulci uu'i'ipiiLCu assuPicu 

1 A« ^ : £ i XX; 

! noi'ipxxtxci scLcxiiya 



Pt = (Et»<30''2 ) )/<2'Zimp ) 
! 



2050 SfrqrPRINT 

28G0 PRINT 

870 PRINT 

380 PRINT 

890 INPUT 



“Input the sampling frequency desired. This frequency 
"determine the range of the sounder. Examples are:” 

SAM. FREQ. = 12500 RANGE = 225 M" 

SAM. FREQ. = 20000 RANGE = 135 M" 

"ENTER THE SAMPLE FREQUENCY" ,Samfreq 



will" 
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2900 
2310 
2320 
2330 
2340 
2350 
23G3 
2373 
2983 
2390 

*7 

3010 
3323 
3033 
3340 
5050 
3050 
3070 
3383 
3093 
3100 
3110 
3123 
3130 
3143 
3153 
3150 
3170 
3180 
3130 
3200 

3210 3peed:INPUT "Enter the tenperature (celalus) “ ,Tenp 



*70 n nk 


Spd»C0kel vi n» ( SQR( 1 +( Tenp/273 ) ) ) 




U iL 


Rspd=Spd/2 




3240 


Rdis t=Tninc*Rsp^ 




3250 


MaArng=*Cnt »Rdist 




3250 


Exrng=( ( INK Maxrng/33 ) )+ . 5 ) *30 . 0 




Toor* 
i '<J 


Far=l . /( Tninc *Spd ) 




3283 


Larobda=Spd/Frcq 




3238 


K3=((2'PI )/Lanbda)Kl/3) 




“7 “700 


Pi 3_lng=( Rspd^Nuro )/Fraq 




“TT 1 O 
1 '<J 


CALL A 1 1 a n Li a 1 1 o n ( F r c q , T a ro p , A t tan ) 




00^3 


! 




j3o3 


! 




3 j43 


PRINT "Enter the relative intensity 


division level. Thi 


3350 


PRINT "value is used to deternine the plot intensity by 


3350 


PRINT "dividing the block average sun by this nunber. “ 


3370 


PRINT "This value is dependent upon 


the gain of" 


3330 


PRINT “the device and will probably 


need adjustpient" 


3330 


PRINT "during run. Start with a value of about 4000" 


3400 


PRINT Pt 





GCLEAR 

Tmnc=l ./Sanfreq 
Tincnt=1000»INT<Tninc/l .E-S ) 
Tine$="TIME ”&yAL$( Timcnt > 
Cnt=UAL(Count$ > 

C3kelvin=331 .5 
Num=iJAL<N$ ) 
Freq=UAL<Freq$)»1000 



I Computation of the D-C offset prior to program run 
! 

! 

Of f_ave= 1 
Of f set=0 
FOR 1=1 TO 13 

CALL Read_ad200< Oat < » ) .Counts .TimeS ,0f f_ave ) 

FOR K = 1 TO (Cnt-1 ) 

Samave=Samave+Oat ( K ) 

NEXT l< 

PRINT “COMPUTING D-C OFFSET" 

NEXT I 
Of f_ave=0 

Of f 3et=Samave/( 10*(Cnt-l > ) 

PRINT "OFFSET IS .Offset 

! 

computation of the speed of sound at a given temp and the 
range of detection of the device 
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3410 

3420 

”7 4 

3440 

j450 

j4S0 

3470 

3480 

3430 

3580 

3510 

3520 

3530 

3540 

3550 

3580 

3570 

3580 

3530 

3S00 

3810 

3820 

3830 

3540 

3650 

3660 

3670 

3680 

3630 

3700 

3710 

3720 

3730 

3748 

3750 

3760 

3770 

3780 

3730 

JO'O'O 

-7n 1 
JO I O 

“7 nn r\ 
JO*, o 

■rn-7r\ 
JO JO 

3040 

3850 

3860 

3870 

3830 

3830 

3300 

3310 



Lvl: INPUT Ilvl 
I 
! 

! set up the plot 



AgainrCALL Plot_3etup(Nplot ,Site$ ,Ma>.rng ,Scale_y .Exrng ,Zone$ ) 
Nrec=0 

FOR 1=1 TO 200 
Ct< I )=0 
NEXT I 
Y=0 

Npoint 1 =0 
OUTPUT KBD;"L"i 
! 

I 

Sync: ! synchronize data collection with clock 
T1=INT(TIME0ATE NOD 86400) 

IF TKT0 THEN T0=T0-86400 
IF T1-T0<1 THEN GOTO Sync 
T0=T1 



Data collection and reduction 



Read_sig: ! read the A-D converter 

CALL Read_ad200<Dat< ♦ ) .Counts ,Tine$ , Of f_ave ) 

Npoint=INT<Tl MOD 3600 MOD Nplot) 

Npoint2=Npoint -Npoint 1 

Nun 1 =Nun- 1 

Nun2=Cnt-Nun 

Dvd= 1 . /< Nun ) 

K-0 

I 

! 

FOR 1=1 TO Nun2 STEP Nun 
R=I*Rdist 
IF KFar THEN R=1 
K=K+1 
r\un_ave=0 
FOR J = I TO I I Nun 1 
X=Dat( J )-0f f set 
Run_ave=Run_ave+X »X 
NEXT J 

Ns=SQR( Run_ave*Dvd ) 

Point<K )=R*(Ns-Corr ) 

IF Kl>10 THEN 

Ct(K )=Ct(l< ) + < (Point(K )''2)»EXP<2*R*Atten) ) 
Ct<K )=<Ct(K)»( (D-SDSE-SI'E )»Er )/Zinp 
END IF 
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3920 

3930 

3340 

3950 

39G0 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

40G0 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

41G0 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

42G0 

4270 

4280 

4290 

4300 

4310 

4 j20 

4330 

4340 

4350 

43G0 

4370 

4380 

4390 



IF Point(K )>ABS(327G7 > THEN Point ( K >='327G7 
D2<K+10)=INT(Point(K)> 

NEXT I 

IF Kl>10 THEN 
Y=Y+1 
END IF 



! Noise correction routine 

I 

! 

IF Ns<(Limit+5) THEN 
Noise=Noi5e t-Ns 
K1=K1+1 
IF Kl>10 THEN 

Corr=Noise/K 1 
Limi t=Corr 

ELSE 

Corr=0 
END IF 
END IF 
! 

I 

I Plotting of the data 



Knax=K 

REDIM D2(Krnax+15) 

! 

! positioning the data on the plot by tine of trace 
! 

Tine=( Npoint2/Np lot )»420 
Tinedist=2+((Npoint2+Npointl)/Nplot)»420 
IF Npointl=0 THEN Tine=5 
IF Npoint<G THEN Tinedist=G 
End_pt=Scale_y »2G0 
D2(2 )=Kna;< 

D2(3 )=INT(Tenp ) 

D2(4 ) = INT(Tine) 

D2(5)=INT(Tinedist»10) 

UINDOU 0,420 ,O,2G0 
GRAPHICS ON 

! set the vertical increnent 

Inc = ( ( Nun*Tninc »Rspd»End_pt )/Maxrng ) 
D2(G)-INT(Inc«100O) 

! 
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4400 


1 compute the intensity of the return, move to 


the prop 


4410 


1 coordinates and plot the appropriate co,lored 


block 


4420 


1 




4430 


FOR K»1 TO Kmax 




4440 


Z=Point(K)/Ilvl 


- 


4450 


IF Z>1 THEN Z=1 




4450 


IF Z<0 THEN Z=0 




4470 


AREA INTENSITY 1,1,1 




4480 


Dis=INT(K»Inc )+l 




4490 


MOOE Timedist.Dia 




4500 


RECTANGLE Time .Inc .FILL 




4510 


NEXT K 




4520 


1 




4530 


! keep an account of the trace numbers taken on the pl< 


4540 


1 




4550 


Nrec=Nrec+ 1 




4550 


Rec_num*INT( Nrec ) 




4570 


D2( 1 )=Rec_num 




4580 


Npoint 1 *Npoint 




4590 


1 




4600 


I 




4610 


! save routine - if function key is set then 




CD 


! the plot will be saved 




4630 


I 




4640 


1 




4650 


IF Cntrl-1 THEN 




4660 


ASSIGN SFilel TO FilelS 




4670 


OUTPUT QFilel ,Rec_num;D2( * ) 




4680 


END IF 




4690 


1 




4700 


i 




4710 


! graphics dump of plot after IS minute intervals 


4720 


! 




4730 


! 




4740 


IF Tioiedist>415 THEN 




4750 


Prt_dmp: PRINTER IS 701 




4760 


PRINT “ 




4770 


PRINT " 




4780 


PRINT " 




4790 


Plotnum=Plotnum-i ' 




A OirKOt 


DUMP GRAPHICS 4701 




4310 


FOR 1=1 TO Kmax 




4 r\n n 
404. It; 


Ct(I )=Ct< I )»<:<Temph273)"2)«( I ./,0038)*K3» 


( 1 . /Y ) 


4830 


Ct( I )=Ct< I)»( 1 ,/( A»G ) ) 




4840 


Ct(I )=Ct< I )»( 1 ./Pt ) 




4850 


Ct(I )=Ct( I )/Pls_lng 




4860 


IF Ct( I XI .0E-'r00 THEN Ct (I )= 1 . 0E- 1 00 




4870 


NEXT I 




4880 


! 
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4830 

4900 

4910 

4920 

4930 

4340 

49S0 

49G0 

4070 

4380 

4090 

5000 

5010 

5020 

5030 

5040 

5050 

5080 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5180 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5280 

5270 

5280 

5200 

5300 

5310 

5320 

5330 

5340 

5350 

5380 

5370 



CT'‘2 Conputat ions and plats 



1 

GCLEAR 

yiEWPORT 15,120,15,80 

WINDOW 0,300 ,0,Exrng 

AXES 30 ,Exrng/ 1 5,0 ,0 ,30 ,Exrng/3 

CLIP OFF 

CSIZE 2, .8 

LORO 8 

FOR M=0 TO 300 STEP 30 
MOOE M,-Exrng/45 
LABEL "1 .0E";(<M/30)-10) 

NEXT M 

MOUE 140,-15 
CSIZE 4 

LABEL “LOG OF CT“2" 

LORG 8 

FOR M=0 TO Exrng STEP INT(Exrng/3) 
MOWE 0,M 
LABEL M 
NEXT M 
LOIR PI/2 
LORG 8 

MOOE -40,Exrng/2 
CSIZE 4 

LABEL "RANGE IN METERS" 

LOIR 0 
LORG 4 

MOOE 150,Exrng+3 

LABEL "TIME AL'ERAGED CT"2 ";Site$ 

CLIP ON 

MOUE 300,0 

FOR 1=1 TO Knax 

Yp 1 = < Maxrng/Kmax ) * I 
Pnt = (LOG(Ct( I ) )/2. 3025851 ) + 10 
IF Pnt>10 THEN Pnt=10 
IF Pnt<0 THEN Pnt=0 
Pnt=Pnt »30 
DRAW Pnt,Ypl 
NEXT I 

DUMP GRAPHICS 5701 
PRINT " 



PRINTER IS CRT 

GCLEAR 

Kl = 10 

Limi t=Corr 
Noise=10»Corr 
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5380 IF Cntrl=l THEN 

5390 ASSIGN SFilel TO » 

5400 END IF 

5410 Cntrl=0 

5420 ! 

5430 ! 

5440 ! if the plot is to be saved, the file is created 

5450 ! and named 

54G0 ! 

5470 ! 

5480 IF Save_plt=l THEN 

5430 CALL File_init(Disc_address$ ,Nrec ,Filel$ ) 

5500 Save_plt=0 

5510 Cntrl=1 

5520 END IF 

5530 ! 

5540 ! 

5550 ! start the next 15 minute plot 

5580 ! 

5570 ! 

5580 GOTO Again 

5590 END IF 
5600 ! 

5610 I 

5620 ! next trace 

5630 ! 

5640 ! 

5650 GOTO Sync 
5660 ! 

5670 ! 

5680 ! toggle set if plot is to be saved 

5690 ! 

5700 ! 



5710 Svpt: Save_plt=l 
5720 GOTO Read_sig 

5730 ! 

5740 ! 



completion routine 



5790 
5800 
5810 
6 000 
6010 
6020 



END 

! 

! 

! SUBROUTINE SECTION 
! 

! 
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G030 SUB Freq_init(Freq$ ,N$ ) 

G043 ! 

6050 ! 

6060 f SETUP OF THE HP3314A FUNCTION GENERATOR 
6070 ! 

6080 ! 

6030 • Ap$»"AP* 

6100 Frq$=“FR" 

6110 Np>$=”NI1" 

6120 En$="EN" 

6130 UoS=>"UO" 

6140 H2$=“KZ* 

El 50 INPUT “FREQUENCY DESIRED <kHZ) (5 RECOMMENDED )“ ,Freq$ 

6160 INPUT “AMPLITUDE DESIRED (U . .1.5U MAX) “ .Amp$ 

6170 IF UAL(Amp$»l .5 THEN 

6180 Amp$=»“1.5“ 

6190 PRINT “AMPLITUDE OF FUNCTION GENERATOR IS 1.5 U“ 

6200 END IF 

6210 INPUT “NUMBER OF CYCLES PER BURST (INTEGER) <100 RECOMMENDED )“ ,N$ 

6220 OUTPUT 707j"MO3“ 

6230 OUTPUT 707;“SR2“ 

6240 OUTPUT 707;Ap$&Amp$&Uo$lFrq$8.Freq$&Hz$8.Nm$8,N$&En$ 

6250 SUBEND 
6260 I 
6270 ! 

6280 ! 

6500 SUB Init_ad200 
6510 ! 

6520 ! 

6530 ! INITIALIZATION OF THE A-D CONVERTER 

6540 ! 

6550 ! 

6560 Ad_sel_code= 1 7 

6570 Dunny—READIOv Ad^5 a 1 c o da , ) 

6530 URITCIO Ad_3al_code ,0 ; 0 

6530 CONTROL Ad_sal_coda ,0; 1 

6600 SUBEND 
6610 ! 

6620 I 
6630 I 
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7000 

7010 

7020 

7030 

7040 

7050 

70G0 

7070 

7080 

7030 

7100 

7110 

7120 

7130 

7140 

7150 

71G0 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

72G0 

7270 

7280 

7230 

7300 

7310 



SUB Set time(Zone$) 



SET THE TIME DATE RECORDER 



PRINT “UHAT TIME REFERENCE ARE YOU USING? INPUT:” 

PRINT ” 1 FOR UNIUERSAL TIME" 

PRINT “ 2 FOR LOCAL TIME* 

PRINT " 3 FOR YOUR OWN CLASSIFICATION" 

INPUT K 
IF K-2 THEN 

Zones-" (LOCAL)" ' 

ELSE 

IF K-3 THEN 

INPUT “WHAT IS YOUR TIME REFERENCE" .ZoneS 

ELSE 

Zones-" (UTC)“ 

END IF 
END IF 

IF TIMEDATE<OATE( *'14 AUG 1984") THEN 

INPUT "ENTER “*DD MMM YYYY" " " .DateS 
INPUT "ENTER " "HR : MIN : SC " “ " .TineS 
SET TIMEDATE DATEC DateS )+TIME( TimeS ) 

PRINT DATE$< TIMEDATE ) .TIMES(nMEOATE) 

Tstart-TIMEDATE 

T0=Tstart MOD 8S400 • 

END IF 

SUBEND 

! 
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7500 

7510 

7520 

7530 

7540 

7550 

75B0 

7570 

7580 

7590 

7B00 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7730 

7800 

7810 

7820 

7830 

7840 

7350 

7860 

7870 



SUB Read_ad200< INTEGER Dat(») BUFFER ,Count$ ,Txne$ , Of f_ave ) 
! 

! 

! INFOTEK A-D ROUTINE SET UP FOR EXTERNAL TRIGGER 



Ad_sel_code=1 7 
I 

! INITIALIZATION OF THE A-0 CONUERTER 
! 

OUTPUT Ad_3el_code: "RESET" 

OUTPUT Ad_sel_codei "INTERNAL" , "COUNT "&Count$ . "HOLDON" 
OUTPUT Ad.sel.coda; "DELAYON-, "SELECT Isl end" ,TiMe$ 
OUTPUT Ad_sel_code 5 "STATUS" 

ENTER Ad_sel_code ;Resp$ 

IF Reaps-" “ THEN 

ASSIGN @Ad200 TO Ad_sel_code ; WORD 
I 

I triggering of the HP3341A 
! 

IF Off_ave-0 THEN 
TRIGGER 707 
END IF 

ASSIGN @Buf TO BUFFER Dat< * ) 

TRANSFER @Ad200 TO ©BufjWAIT 
OUTPUT Ad_sel_code;“ “ 

OUTPUT Ad_sel_code; "STATUS" 

ENTER Ad_5el_code iRespS 

IF RespSO" THEN ‘ 

PRINT ‘ERROR- “;Re3p$ 

END IF 

ELSE 

PRINT “ERROR DURING INITIALIZATION -"jReapS 
END IF 
SUDEND 
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8QQ0 

8010 

8020 

8030 

8040 

8050 

8060 

3070 

8080 

8090 

8100 

8110 

8120 

8130 

8140 

8150 

8160 

8170 

8180 

8180 

8200 

8210 

8ZZ0 

8230 

8240 

8250 

8260 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8360 

8370 

8330 

8380 

0 A Hi Hi 

8410 

8420 

8430 

8440 

8450 

8460 

8470 

8480 



SUB Plot_setup< Nplot ,Site$ .Maxrng ,Scale_y ,Exrng ,Zone$ ) 
! 

! 

! SET-UP OF THE TIME PLOT OM THE CRT 



Scala_y=Maxrng/Exrng 
GRAPHICS ON 
LINE TYPE 1 
UIEUPORT 15,120,15,80 
WINDOW 0 , Nplot ,0 ,Exrng 

AXES Nplot/15,Exrng/15,0,0,Nplot/3,Exrng/3 
CLIP OFF 
CSIZE 4, .6 
LORG 6 

Tl-TIMEDATE MOD 86400 
Hrs=Tl DIU 3600 
T2=T1 MOD 3600 
Min-T2 DIU 60 
Qtrhr'Min DIU IS 

FOR M=0 TO Nplot STEP rNT<Nplot/3) 

MOUE M,-Exrng/45 
Qtrtnin=Qtrhr» 1 5+( M»3/Np lot )»S 
IF Qtrnin=G0 THEN 
Qtrnin“0 
Hrs=Hr5+ 1 
END IF 

LABEL USING ‘DO ,A ,ZZ* iHra 5 " : * 5 Qtrnin 
NEXT M 

MOUE Nplot/2,-15 
LABEL "TIME "iZone$ 

I 

I LABEL ORDINATE 
I 

LORG 0 

FOR M=0 TO Exrng STEP INT(Exrng/3) 

MOUE 0,M 
LABEL M 
NEXT M 
LDIR PI/2 
LORG G 

MOUE -Nplot/7 ,Exrng/Z 
LABEL "RANGE IN METERS” 

TITLE 

LDIR 0 
LORG 4 

MOUE Nplot/2 ,Exrng+3 
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8490 

8900 

8510 

8520 

8530 

8540 

8000 

9010 

9020 

9030 

9040 

9050 

9080 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9180 
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9200 
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9220 
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9280 

9290 

9300 
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3 -rn 

UtL'O 

G33G 

934Q 

3350 

9350 

9370 

3380 

5390 

9400 

9410 



LABEL DATE$(TIMEDATEr; 
CLIP ON 
SUBEND 



! 






SUB Attenuat ion( Freq ,Tenp ,Atten ) 



! This subprogram calculates the attenuation of the 
! sound in air based upon equations in Neff 1975 
! (source of subroutine: Thesis of R. Fuller) 

I 

! Variables 



I Atom _pres - input atmospheric pressure in mb 
! 

I Atten ~ attenuation coefficient of acoustic wave 



! Att_max - Variable in program. It is the attenuation 
! at the frequency of the maximum attenuation 

I for the given input conditions, 

I 

! F - Ratio of the frequency to frequency at maximum 
I attenuation, 

I 

! Fmax - Frequency of the maximum attenuation 
I 

! H - variable used in the integration of excess attenuation 



Pstar - variable used in intermediate calculations 



I Tstar - variable used as an intermediate in calculation 
1 of attenuation, 

I 

I Uat _pres - Atmospheric water pressure in mb. 

! 

INPUT “Enter the atmospheric pressure in mb" ,Atom_pres 
INPUT "Enter the atmospheric water pressure in mb" ,Uat_pre3 
H=1 0O»Wat_pres/Atom_pre5 
Tstar-( 1 ,8»Temp + 492 )/513 
Pstar=Atom_pres/ 1 01 4 

Fmax-( 10+BG00*H+44400»H»H)»Pstar/T3tar".8 
Att_max».0078»Fmax»Tstar''(-Z.5)*EXP(7,77*( 1-1/Tstar ) ) 
F=Freq/Fmax 

Atten-( Att.max/304.8)»( C. [ 8»F )''2+( 2*F*F/( 1 +F»F ) )"2 )" .5 

Atten»( At ten-M . 74E- 1 0»Freq*Freq )/4.35 

SUBEND 
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9420 

9430 

9500 

9510 

9520 

9530 

9540 

9550 

95B0 

9570 
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9510 

9520 



! 

! 

SUB Fi le_init ( Di5c_address$ ,Nrec .Filet $ ) 

1 

I 

! CREATE THE STORAGE FILE ON THE DISC FOR 
! THE REDUCED DATA 
I 
! 

INPUT "ENTER THE REDUCED DATA OUTPUT FILENAME 
Filet $=File$&Disc_address$ 

CREATE BOAT Fi le t $ . t 80 ,400 ’ 

ASSIGN QFilel TO FileiS 

! 

SUBEND 



.Files 
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APPEHDIX B 



SPEAKER AND ARRAY ANALYSIS 

In the design of our echosounder, it was determined 
that a rapid decay time was required to obtain accurate 
short range information. Weight restrictions involved with 
equipment transportation plus maximum response at high 
frequencies led to our decision to use piezo ceramic 
speakers. The Motorola KSN 1005A speaker was selected 
based on these requirements and the specifications charted 
in the Motorola catalog [Ref. 16] and reproduced below in 
Figures 23 through 25. 



Nominal Power, Impedance, and Distortion Ratings 




Fig. 23. Speaker Ratings 
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Typical Frequency Response 




500 1000 2000 5000 10000 20000 40000 

Fig. 24. Speaker Frequency Response 



Dimensions: KSN 1005A, KSN 1003A 



FOUR 218- 15 5mfT>) DIAMETER HOLES EQUALLY 
SPACED ON A 3 94" (lOO.Imml DIAMETER 8 C. 




61 



The frequency response charts indicate that a maximum 
response for our speakers occurs at a resonant frequency of 
5000 Hertz. This frequency was used as the baseline from 
which all our measurements are made. 

Using the speakers in an anechoic chamber, the average 
e"i voltage decay time was measured to be approximately 900 
usee (Fig. 26). This decay time translates into a sound 
propagation distance of just over 15.0 centimeters (at STP) 
from the speakers. Considering our requirements, this 
speaker is ideally suited to serve our purpose. 




Fig. 26. Speaker Decay Time Trace 
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B. ACOUSTIC ARRAY 



Our next consideration was the echosounder array 
pattern. Ideally, the acoustic sources should bQ placed 
exactly one half wavelength apart. At a frequency of 5000 
Hertz, this would require spacings of 3.4 centimeters (at 
STP) which is physically impossible for the speakers we 
have chosen. The closest possible spacing is 7.62 
centimeters between sources after shaving off the flange of 
the horn (Fig. 25). 

From Kinsler, et al . [Ref. 14], the equation for the 
directionality factor of a simple line array is derived as: 



H(N,6) 



1 

N 



N 

sin ( kd sin 0 ) 

2 



1 

sin ( kd sin 0 ) 

2 



where 

k is the wave number (2n/ A ), 
d is the distance between sources, 

N is the number of sources, and 

0 is the angle measured from a line perpendicular to 
the array to the direction of interest. 

However, this equation assumes simple point sources which 
does not adequately describe our speakers. It was necessary 
to couple this equation to the directional factor for a 
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piston source which is also identified in Kinsler, et al . 
[Ref. 14: p. 108] as: 

I 2 Ji ( ka sin ( 9 ) ) 

D(0) = 

I ka sin ( 9 ) 

where 

k is still the wave number, 

a is the radius of the piston source, 

9 is the angle measured from a line perpendicular to 
the array to the direction of interest, and 

Ji ( X ) is a Bessel function of Order 1 with an 
argument x . 

These two equations were combined together to form the 
equation of directionality for a linear array of piston 
sources, L(N,9), by simple multiplication. 

L(N,9) = H(N,9) • D(9) 

This equation was incorporated into a computer program, 
which may be found in Reference 3, and used to generate 
linear array patterns for various numbers of piston sources 
(Figs. 27 through 31). 

It was concluded from these plots that five linear 
elements would provide the best combination of forward 
directionality, sidelobe suppression, physical size and a 
relatively low cost. Then, in order to enhance both array 
efficiency and symmetry, we settled on a five by five 
element array design with 7.26 centimeter spacing between 
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Fig. 27. Three Element Array Fig. 28. Four Element Array 




Fig. 29. Five Element Array 
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speakers in both the vertical and horizontal directions. 

After verifying the manufacturers polarity designation 
for 35 speakers, we obtained Lissajous plots for each 
individual speaker (Fig. 32). Based on the speaker’s output 
to input voltage ratios as illustrated by these plots, we 
were able to rank all our speakers by signal efficiency. 
It was based upon this criteria that we selected the 25 
most efficient speakers for the array, placing the best 
speakers at the center and subsequent ranking speakers 
further toward the sides and corners. 




Fig. 32. Lissajous Plot 
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The 25 selected speakers were mounted in a five by five 
planar array on a balsa wood insulated bilayered sheet 
metal board. After wiring all the speakers in parallel, we 
surrounded all the electrical connections and speaker backs 
with two 3.0 centimeter layers of foam insulation 
sandwiching a 1.0 millimeter lead sheet. Then the entire 
array mounting was enclosed in a 44 by 44 by 5 centimeter 
sheet metal box. This design was chosen to suppress 
virtually all acoustic energy propagating out the rear 
hemisphere of the array, while shielding the array from any 
external electrical interference (Figs. 33 and 34). 
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Fig. 33. Array Photo 
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Fig. 34. Array Photo 
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APPEMDI3L-C 



ENCLOSURE DESIGN 

Acoustic echosounding has proven to be an extremely 
useful technique for probing and analyzing the lower 
atmosphere. In order to most efficiently utilize the 
acoustic waves transmitted and later received by this 
remote sensing method, it is essential to have an efficient 
antenna with highly directive beams and strongly suppressed 
sidelobes. Antenna design becomes increasingly more 
important in a noisy environment where noise pollution 
within the sidelobes may dominate the desired signal within 
the main lobe. Hall and Wescott [Ref. 12] showed that 
sidelobe suppression improved with higher frequencies. 
Their studies showed that the measured 90 degree sidelobe 
suppression ranged from 38 dB at 1 KHz to 50 dB at 5 KHz. 
Furthermore, any significant improvement in sidelobe 
suppression could only be obtained by surrounding the 
antenna with an acoustic energy absorbing cuff or shroud. 

In an effort to maximize our antenna main lobe to 
sidelobe power ratio, we intend to operate only at high 
frequencies as discussed in the previous section. 
Additionally, we have designed an acoustic energy absorbing 
enclosure. Many designs were considered based upon 
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previous research in the field of echosounding [Refs. 17 
through 20], In addition, we obtained the actual acoustic 
beam patterns for our array using a computer program 
written by LCDR Butler [Ref. 21] which we modified for our 
purposes. This modified version of LCDR Butler’s program 
may be found in Reference 3. By rotating the array in an 
anechoic chamber, we were able to produce highly accurate 
polar plots of the array beam patterns (Figs. 35 and 36). 



5X5 ACOUSTIC ARRAY BEAM PATTERN 




R=3.60 M MIKE VT= 8.2465 INPUT-5.0 V FREQ-5.0 KHZ 
Fig. 35. Polar Plot of 5 X 5 Array Beam Pattern 
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DIAGONAL 5X5 ACOUSTIC ARRAY BEAM PATTERN 




R-3.60 M MIKE VT- 3.9497 INPUT-2. S V FREQ-5.0 KHZ 



Fig. 36. Polar Plot of 5 X 5 Array Beam Pattern 

at a 45® Aspect 
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These polar plots and the actual coordinates 
corresponding to the individual data points confirmed the 
computer prediction for the five element line array (Fig. 
29). Additional polar plots obtained by varying both the 
input array voltage and range between array and microphone 
further support our claim that for a 5 KH 2 carrier 
frequency, the main lobe is confined to a divergence angle 
of 20 degrees. Since it is our aim to suppress all 
sidelobes and utilize solely the main lobe, we chose not to 
taper our enclosure as most previous researchers had. 
Rather we designed the enclosure based upon the dimensions 
of the array itself and the acoustic beams it generated. 

Plywood was used for the construction of the enclosure 
and provided not only a rigid, inexpensive framework, but 
also proved to greatly attenuate external noise 
interference. Anticipating all kinds of weather conditions 
during data collection, the plywood enclosure was first 
waterproofed with four coats of marine varnish. Grooved 
joints, caulking and weather stripping were also design 
considerations . 

A millimeter layer of lead can suppress an acoustic 
signal as much as 40 dB (Figs. 37 and 38). About 7.0 
centimeters of corrugated, egg-carton design foam can 
suppress a signal another 3 to 4 dB (Figs. 37 and 39). 
Together they make an extremely efficient absorbing 
material for use in our enclosure. 
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Fig. 



37. Response Reference, No Insulation 



START 4 00 a Hz RAMGE: -51 dBV STATUS: PAUSED 

AttlAG RHS:25 




X: 5000 HZ Y;-60.77 dBV 



Fig. 38. Signal Suppression by Lead 
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Fig. 39. Signal Suppression by Foam 



Lead-lined absorbing foam was glued to all inner 
surfaces of the enclosure with two 1.0 millimeter layers of 
lead overlapping at all corners. Strong aluminum brackets 
were used to connect the four side panels to each other as 
well as to the enclosure base (Figs. 40 through 42). 
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Fig. 40. 



Enclosure Photo, Fully Assembled 
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Fig. 41. Enclosure Photo, Interior 
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Fig. 



42. Enclosure Photo, 



Base with Array 
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APPENDIX D 



ECHOSOUNDER OUTPUT 

Sixteen echosounder output traces were included to 
exhibit typical atmospheric activity. Many of the plots, 
such as the 14:45, 15:30, and the 17:30 of 26 April and the 
10:00, 11:15, 13:00, 14:15, 15:15, 16:15, and the 17:15 of 
27 April, have convective plumes which are prevalent 
whenever a heat flux between the surface and atmosphere 
exists. The plots of 18:15 and 19:15 on the 26th of April 
and of 18:00, 18:45 and 20:00 on the 27th of April clearly 
show the passage of the neutral event which is encountered 
when the atmospheric and surface temperature difference 
becomes negligible. Finally, the plot of 16:30 on 26 April 
can be associated with strong winds which exhibit a 
somewhat uniform return for all altitudes across the entire 
15 minute interval of the trace. 

The horizontal axis labels of the Ct 2 plots are small 
and may be hard to read. This axis is a log scale 
beginning with l»10-io at the left side of the plot and 
ending with 1*100 at the right side of the plot. Each tick 
mark moving left to right along this axis represents an 
integer increase in the exponent of 10 . 
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RRriGE IN METERS RRNGE IN METERS 
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TIME (LOCAL) 



TIME AVERAGED CT^2 NPS ROOF 
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RRt-JGE irj METERS RRNGE IM METERS 



26 Rp r 1987 NP3 ROOF 




TIME (LOCRL) 



TIME AVERAGED CT-2 NFS ROOF 




LOG OF CT'^Z 
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RRfIGE IN METERS RBNC.E IN METERS 



2b Apr 1987 NFS ROOF 




16:30 18:35 18:40 18:45 

TIME (LOCAL) 



TIME AVERAGED CT-2 NFS ROOF 
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